Competing effects of interactions and spin-orbit coupling in a quantum wire 
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We study the interplay of electron-electron interactions and Rashba spin-orbit coupling in one- 
dimensional ballistic wires. Using the renormalization group approach we construct the phase dia- 
gram in terms of Rashba coupling, Tomonaga-Luttinger stiffness and backward scattering strength. 
We identify the parameter regimes with a dynamically generated spin gap and show where the 
Luttinger liquid prevails. We also discuss the consequences for the operation of the Datta-Das 
transistor. 

PACS numbers: 71.70.Ej 
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The main goal of recent studies in the field of spin- 
tronics is to invent ways of manipulating the electron 
spin with an efficiency comparable to that of present- 
day electronics (which manipulates charge) Q. One of 
the elementary spintronic devices, the Datta-Das tran- 
sistor 0, has been proposed more than a decade ago, 
its basic ingredient being a ballistic quantum wire with 
sufficiently strong Rashba spin-orbit interaction Q . The 
latter is required for creating a sizeable spin precession. 
Depending on spin orientations in the source and in the 
drain one can modulate the current flowing through the 
device and thus implement in principle ON/OFF states. 
The strength of the spin-orbit couplingcan be tuned by 
applying a gate voltage to the system 

To understand the feasibility of such a device as well 
as its basic operation it is important to investigate the 
effects of spin-orbit coupling on both transport and mag- 
netic properties of (essentially one-dimensional) interac- 
ting electrons. Recently some progress towards a so- 
lution of this problem has been made on the basis of 
the Tomonaga-Luttinger model 0. In the present ar- 
ticle we show that the a priori assumption of validity 
of this approximation is not always justified. We find 
in fact that the combined effects of Coulomb interaction 
and Rashba coupling can generate a spin gap and thus 
radically change the physical characteristics of the Datta- 
Das transistor. In particular, we establish the parame- 
ter range where the correlation functions decay exponen- 
tially, and the Datta-Das device becomes non-operating. 

We consider a narrow ballistic wire described by a one- 
dimensional model of interacting electrons, including the 
Rashba-type spin-orbit coupling, Vr = an<7 y p x , where 
an is the coupling strength and o~ y is a Pauli matrix. 

The electron-electron interaction (assumed to be weak) 
may be decomposed into the different scattering pro- 
cesses 0. We include not only the standard forward 
scatterings of the Tomonaga-Luttinger model (coupling 
parameters gi\\, g2\\, 92±, <74||, <?4_i_), but also the backward 
scattering term (gi±). At the same time we do neglect 
Umklapp processes (gs) since the system considered here 
is far away from half-filling. 



Although in the absence of spin-orbit coupling (and 
for repulsive electron-electron interactions) the backscat- 
tering term is usually irrelevant, this is not always the 
case in the presence of spin-orbit coupling. As will be 
shown below, the combined effects of these two processes 
together with the forward scattering lead to the dynam- 
ical generation of a spin gap in the excitation spectrum 
in a wide range of coupling constants. In this range 
the ground state phase diagram qualitatively differs from 
that of the Tomonaga-Luttinger model used before 0,0). 

We assume weak bare couplings and therefore lin- 
earize the single-particle spectrum around the two Fermi 
points ±kp. Decomposing the field operators in the 
standard way 6] into right (r = +) and left movers 
(r = — ), ip„{x) = J2 r ^r,(r(x) exp(irk F x), we obtain 
the continuum limit of the fermionic Hamiltonian H = 
Hq + Hr + Hi n t, where Ho and H mt are the usual one- 
and two-particle parts, respectively, while the spin-orbit 
term reads 



H T 



-ia R k F ^rj cb:(^(x)^(x)-^l >L (x)'ip r ^(x)). (1) 

We remark that terms involving gradients of the fields 
ip rir7 (x) have been neglected because they would gen- 
erate irrelevant operators for non-vanishing electron- 
electron interactions ||. All parts of the Hamilto- 
nian can be bosonized in terms of fields 4> c {x),9 c {x) 
for the charge and <p s (x),8 s (x) for the spin degrees of 
freedom, with commutation relations [4>\(x),9\i(x')] — 
(—i/2)5\\>sigii(x — x') 0. The Rashba term 



H R = 



-J sin(v^2~7T0 s ) sin(v / 2~7r0 s ) 



(2) 



(with a short distance cut-off a) involves only <p s and 
6 S and therefore does not affect spin-charge separation 
(the same holds for the backscattering term). The ki- 
netic energy and the forward scattering terms become 
H = H c xi + H s q, where the charge part H c q and the 
free spin part i? s ,o both have the familiar Tomonaga- 
Luttinger form, 



Hx, = y J dx [K^Jx? + -^{d x <Px?\ 



(3) 
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The parameters u\, K\, A = c, s depend in a simple way 
on the coupling constants gi,\\, 92,\\, 92,±, 9i,\\, 9i,± (EH- 
Treating H R as a perturbation we notice that its vac- 
uum expectation value vanishes (a consequence of non- 
zero conformal spin) . Thus the lowest-order effect of the 
Rashba one-particle process is seen to be absent, even 
though the process is strongly relevant. One then has 
to study the higher-order expansions in order to find 
contributions with zero conformal spin (cf. Ref. 0,E3) 
The second-order contribution (proportional to a R ) cor- 
responds to effective two-particle interactions 



(4) 



a first term describing a backscattering process already 
included in the Hamiltonian and a second term repre- 
senting a spin-nonconserving process which has not been 
met before. We therefore have to add the new interaction 
term in the renormalization procedure, with a vanishing 
initial coupling constant. In this way the operator prod- 
uct expansion governing the renormalization group (RG) 
flow is closed. We also notice that the process in ques- 
tion usually emerges in the field-theoretical description 
of systems with completely broken spin-rotational sym- 
metry |lllll4j. 

The bosonized Hamiltonian for the spin part can be 
written as H s = H St o+H R +Hbf, where H s q and H R have 
been given before and -f/bf includes the backscattering 
and spin flip parts, cf. Eq. 10} , 



bf : 
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1_L 



(27ra) ; 



- / dx cos(V8ft(f>s) + 7—^~To f dxcos(V8n9 s ). 
J (2Trayj 

(5) 



The Hamiltonian H s has the dual symmetry 



0„, K„ 



k: 



9ix 



9f, 



(6) 



which will be used below to obtain the phase diagram 
of the model in the full parameter range. Thus, if a 
transition occurs for some set of parameters there must 
be also a transition for the dual set of parameters. 

Standard perturbative RG analysis 12] up to the third 
order in couplings yields the set of equations 

= {4-{K. + K^ 1 ))y R + y ± y f y a , 

^ = 2(l-K s ) y± + (K- 1 -K s )y R +^(yl-yj), 

^ = 2(1- K-')y f -(K-i-K s )y R - y -l (yl-y% 



dK s 
dl 



5 (tf 



(7) 



where I measures the logarithm of the length scale, 
V± = gi±/mt s , y f = gf/nu s and y R = (2a R k F a/iru s ) 2 . 
Note that it is y R that enters J7J|, not a R . The RG 
equations are therefore independent of sign(a^), a con- 
sequence of time-reversal symmetry. The initial condi- 
tions are K S (Q) = K , y±(0) = y± <0 , y R {0) = y R ,o and 



y/(0) = 0. Similar RG equations have been derived in the 
context of spinless fermions on a ladder HEii. The RG 
equations Q have three weak-coupling fixed points: (1) 
y* ± = yj = y* R = and K* arbitrary ; (IIa,b) A"* = f, 
V*R = Vj = ^Vx an d Vx arbitrary. The spin- wave fixed 
point 13] (I) corresponds to the noninteracting system 
with no spin-orbit coupling. The fixed points (II a,b) cor- 
respond to the critical Ashkin- Teller (AT) model fi a ll ^ ]. 

A rough idea of the RG flow can be given by assum- 
ing the spin stiffness K s to be constant {=Kq) and by 
neglecting cubic terms in J7J). The solutions are 



Vr(1) 

yx,(f)(i) 



VR,oe 

l -L,(/) 



l7R 



A 2l[l-K 



Bx,( f) e 1 ^, 



(8) 



where B ±t(f) = ±y Rfi [K^ - K ]/[2 ± A T A^ 1 ] and 



A 



x,(f) = yx,(f),o - -B_l,(/) 



The one-particle Rashba pro- 



cess is relevant if the exponent 7^ = [4 — (A + K^ 1 )] 
is positive, i.e. for 2 — V3 < Kq < 2 + \/3. The 
amplitude y R ,o is larger than both B± and Bf for 
(V5-l)/2 < A < (v / 5+l)/2. This (Rashba dominated) 
region will be investigated more accurately below. For 
A'o < (s/5 — l)/2, backscattering (y±) dominates, while 
for A > (VE+l)/2, spin flip processes (yf) prevail, in 
agreement with the duality relations iffijl. 

We return now to the full RG equations 0. Numer- 
ical solutions for various initial conditions in the regime 
of dominant spin-orbit coupling are illustrated in Fig. ^ 
as flows in the (y R , ^lnA s ) plane. The results confirm 
that K s remains essentially constant for yj_.o = and is 
weakly renormalized for y_i_ t o — 0.4. On the other hand, 
the Rashba coupling y R rapidly grows from its initial 
value y R fi up to some value of order 1 at a length l R , while 
the couplings y± and yf remain small. It is then possible 
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FIG. 1: Flow diagram for the (j/H,|lnA 3 ) plane for yn t o = 
0.07 and various initial values Kq. Full lines correspond to 
y± t o = 0, dashed lines to j/j_,o = 0.4. 

to add a second step in the RG procedure by treat- 
ing H R at l R non-perturbatively through the canonical 
transformation Vv,± = 771 (Vv.TJ. — *V-V.it)- The trans- 
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formed single-particle Hamiltonian is 

H + H R = J dx iply(x)(rv F d x + r'atRk F )ip r y(x). 

Its spectrum consists of two subbands split horizon- 
tally by 2ajik F /v F , where v F = u s K s is the renor- 
malized Fermi velocity. The canonical transformation 
leads to the same type of interaction terms as in the 
original Hamiltonian, with coupling constants yx,f ~ 
K,-l + (y ±J - y/,_i_)/2, K s - 1 « (y ± + y f )/2 (for 
\K S — 1| -C 1; otherwise the relations are slightly more 
complicated). 

Bosonization can be applied to the new Hamiltonian 
in terms of fields (f> s and 9 S . The Rashba term leads to 
a contribution (y/2/n)aRk F J dxd x 9 s which can be ab- 
sorbed in the single-particle part by the shift 



Q s (x) = 9 s (x) + k R x, ®„(x) = (/> a (x), 



(9) 



where kR — \f2a.R.k F /(iru s K s ). The spin part of the new 
Hamiltonian reads 



/ dx\K s {d x Q s Y + -~r{d x ®s) 2 



(10) 



cos(v / 8^$ ;j ) - cos(V8^e s - k R x) 



Vf 



7ror 



For K s > 1 the ordering in the S sector and the Rashba 
process are competing, but since for yR « 1 and yf <C 
1 the last term is wiped out by strong oscillations, the 
subsequent flow involves only the parameters K s and y± . 
The equations of the second RG step have the form of the 
Kosterlitz-Thouless flow 



dK s 
dl 



1 



dy± 
dl 



2(1- K s )y ± , 



(11) 



with the initial conditions K s (0) = l+^(y±(lR)+yf(lR)), 
y±{0) = K.(l R )-l + ±(y ± (l R )-y f {lR)). Eqs. dl imply 
two distinct scenarios, for |yj_| > 2[K S (0) — 1] a flow to 
strong coupling and for \y±\ < 2[K S (0) — 1] a flow to 
the weak-coupling fixed line y*^ — 0, K* (non- universal 
renormalized value). Clearly the line y± = yf, K s — 1 
is critical for arbitrary yR. We find systematically the 
weak-coupling case (for y±^ > 0). Thus, if the flow is 
dominated by the Rashba term, the system is a Luttinger 
liquid, with enhanced spin precession. 

We discuss now the parameter regions where the 
Rashba term is not dominant, and the RG flow is gov- 
erned by Eqs. J7J. We find numerically that for K s ^ 1 
either y± or yf are renormalized towards strong coupling, 
y± for K s < 1, yf for K s > 1, in agreement with the 
approximate analytical solutions ©■ Both cases scale 
to the strong-coupling regime of the sine-Gordon model 
and therefore must have a spin gap A s in the excitation 
spectrum. At the same time, the Rashba term and the 



band splitting are dynamically suppressed. This is an un- 
expected new result, because in the absence of the bare 
spin-orbit coupling there is no spin gap for y± > 0. 

Our findings can be interpreted in terms of two 
commensurate-incommensurate transitions |l5j . In the 
incommensurate phase the Rashba term dominates and 
the effective field theory is equivalent to a Tomonaga- 
Luttinger model with non-universal spin stiffness and 
momentum shift (|5J). In the commensurate phases the 
spin excitations have a gap, produced by backscattering 
and spin flip processes, respectively. Fig. shows the 



C-IC 




0.5 J R,0 



C-IC 



FIG. 2: Phase diagram in the parameter space of initial values 
j/h,o and 2 hi(-Ko) for y± t o = (full lines), y± t0 = 0.4 (dashed 
lines), and for y/,o=0. The figure exhibits clearly the dual 
symmetry of Eq. © . Spin gaps exist both in the phase with 
dominant superconducting correlations (SS or TS(x)) and in 
the spin (charge)-density wave phases (SDW(x) or CDW). In 
the Luttinger-liquid phase (TL(inc)) all correlation functions 
decay algebraically. The TL phase is separated from the two 
spin-gapped phases by commensurate-incommensurate tran- 
sitions (C-IC lines). The bold dotted line (AT at K s = 1) 
corresponds to the critical Ashkin- Teller model. 

phase diagram in the parameter space (yR,o, \ hi-Ko) for 
the initial values yx,o = 0, 0.4 and yf$ = 0. The straight 
line IyiKq — (fixed points (II)) represents the self-dual 
line for yj_,o = ?//,o = 0. The region of the incommensu- 
rate Luttinger-liquid phase widens as a function of yR$. 
Outside of this region, the mean values (<^> s ) or (9 S ) are fi- 
nite, the former for Kq < 1, the latter for Kq > 1. In the 
4> s -ordered phase the dominant correlations are singlet 
superconductivity (SS, for K c > 1) and charge-density 
waves (CDW, for K c < 1). In the ^-ordered phase 
the dominant correlations are the x-component of spin- 
density waves (SDWz, for K c < 1) and the x-component 
of triplet superconductivity (TS^, for K c > 1). The tran- 
sition lines in Fig.[2havc been determined approximately 
on the basis of the RG flow. We note that the exact 
characterization of these commensurate-incommensurate 
transitions would require a non-perturbative analysis, 
which goes beyond the present perturbative RG scheme. 
Spin precession is described by the correlation function 
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f{x) = i<(^ T (.T)+^(x))(^(0) + Vi(0))). (12) 

For a narrow, ballistic quantum wire connecting a source 
at x = to a drain at x = L the quantity |/(L)| 2 mea- 
sures the probability for a particle entering the drain to 
have the same spin orientation as one leaving the source. 
Within the Tomonaga-Luttinger model f(x) is found 
to vary as where 7 = (1/A)(K C +1/K c +K s + 1/K a ) 

depends on the electron-electron interaction through the 
charge and spin stiffnesses K c and K s , respectively. The 
same behavior is expected to occur in our case within 
the Luttinger-liquid phase except that K s is replaced by 
the (non-universal) fixed-point value K*. In principle, 
the presence of irrelevant couplings yj_ and yf can lead 
to extra multiplicative corrections [6j to fix), but they 
appear to vanish identically in our case. In the spin- 
gapped regime where the ordering of the </> s field implies 
the disordering of 8 S field for K s < 1 (and vice versa for 
K s > 1), the function f(x) is expected to decay as the 
corresponding correlator in the sine-Gordon model (lif . 
namely fix) ~ :r _ ( A+1 ) exp(— x/£ s ), where £ s = hvp/A s 
is the correlation length, and A = i\fK~ c — l/v^c) 2 /^- 

We note that the spin stiffness K s is an important 
parameter in the problem. In 5J7(2)-invariant models 
with repulsive interactions K a scales to the fixed-point 
value K* = 1. This is no longer true for systems with 
spin-orbit coupling where the spin-rotation symmetry is 
reduced to t/(l). 

As an example we consider a narrow InAs quantum 
wire with to* = 0.023m e and a R = (0.6 - 4) x 1CT 11 
eVm We choose the wire width d — 5 nm and the 
Fermi wave vector kp — 0.5 x 10 8 m _1 . In this case the 
assumption of a single occupied subband is justified. The 
parameter yn^ depends on the cut-off length a, for which 
a natural choice is the width d. Using these values and 
u s w vf we find that yn,o ranges from 10~ 4 to 10~ 2 . 
In order to estimate Kq, we use the well-known relation 
between the spin stiffness and the Fourier transform 
V e f f (q) of the effective interaction potential for quantum 
wires. The latter is taken in the form proposed in Ref. 
[13 . Thus we obtain ±lnK « 0.15. According to Fig.H 
this corresponds to the spin-gapped phase. The standard 
procedure for the sine-Gordon model yields a value of 
A s in the range (0.01 — 0.1)ef, where ep ~ 4 meV for 
given to* and kp- For temperatures T < A s /fcs = 0.5 — 5 
K and wire lengths exceeding the correlation length, L > 
£ s = (0.4 — 4)/im the phenomena described above will 
play an important role. These regions can be reached 
in present-day devices (L » 2 — 6/im), and it should in 
principle be possible to detect signatures of the spin gap. 

If the material parameters can be tuned close to the 
commensurate-incommensurate transition from a Lut- 
tinger liquid to a phase with a spin gap, the dependence 
of the Rashba coupling on the electric field strength may 
allow to drive the system from one side of the transition 



to the other. Such a control of a spin gap by a gate volt- 
age would represent a spectacular novel field effect. The 
observation of such a subtle phenomenon would of course 
not only be of fundamental interest, but it could also pave 
the way for the fabrication of new types of devices. 

In conclusion, we have found that the interplay be- 
tween electron-electron interactions and Rashba spin- 
orbit coupling in a narrow wire generates a spin gap A s 
for a certain range of parameters. This restricts the op- 
eration of the Datta-Das transistor. In the Luttinger- 
liquid phase, where correlation functions fall off accord- 
ing to power laws, the device may work, but in the spin- 
gapped phase the spatial coherence of spin precession is 
suppressed exponentially, and the device efficiency tends 
to zero at lengths greater than the correlation length £ s . 
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